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Abstract. A simple thermal model with single freeze-out and flow is used to analyze the ratios 
of hadron yields and the hadron transverse-mass spectra measured in y/J^ =130 GeV AuH-Au 
collisions at RHIC. An overall very good agreement between the model predictions and the data is 
achieved for all measured hadron species including hyperons. 



INTRODUCTION 

The main features of the soft hadron production at RHIC, such as the ratios of hadron 
abundances, the transverse-mass spectra, the elliptic flow, or the HBT radii may be 
efficiently understood in the framework of a statistical model which combines the 
standard thermal analysis of the hadron ratios with a suitably parameterized expansion 
of matter at freeze-out [1, 2, 3]. In the present paper we outline the main assumptions of 
the model and concentrate on the discussion of the particle ratios [4, 5] and the hadron 
transverse-mass spectra [6, 7, 8, 9]. The subsequent paper in these Proceedings [10] 
contains the analysis of the elliptic flow and the pion HBT correlation radii. 

The main ingredients of our approach are as follows: i) the chemical freeze-out 
and the thermal (kinetic) freeze-out occur simultaneously, ii) all hadronic resonances 
are included in the calculation of both the hadron yields and the spectra, and iii) the 
freeze-out hypersurface and the flow at freeze-out are defined by the simple expressions 
inspired by the Bjorken model [11]. Below we discuss in more detail these three points. 



Freeze-out 

Our approach includes a complete set of hadronic resonances in both the calculation 
of the hadron abundances and the calculation of the hadron spectra. Then, it turns out 
that the distinction between the two freeze-outs [12] is not necessary. Our analysis 
showed [4] that the decays of the resonances, which are initially present in a heat bath 
at the temperature of 165 MeV, effectively lower the inverse slope parameters of the 
spectra by about 30-40 MeV. This is just the desired effect which explains the typical 
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difference between T^^^^ and T^^^^ (i.e., between the temperature at the chemical freeze- 
out and the thermal freeze-out, respectively). As a consequence, we have found, at least 
for the RHIC data, that no extra elastic rescattering is required in order to describe 
simultaneously both the ratios and the spectra. In other words, we assume one universal 
freeze-out taking place at [1] 

^chem — ^kin = (1) 

Recently, experimental hints have been found in favor of our assumption (1). A success- 
ful reconstruction of the K* (892)*^ states by the STAR experiment [13], together with the 
very good agreement between the measured yield of K*{S92)^ and the prediction of the 
thermal model [3, 14] suggests a picture with the short expansion time between the two 
freeze-outs. Such a picture is natural if the production of particles (hadronization) occurs 
in such conditions that neither elastic or inelastic processes are effective. An example 
here is provided by the sudden hadronization model of Ref. [15]. 



Resonances 

For thermal systems the contribution from high-lying (heavy) states is damped by 
the exponential factor, exp(— m^/T). This fact, at the first sight, suggests that most of 
the resonances present in a hadron gas at a moderate temperature may be neglected. 
However, although the high-lying states are suppressed, their number increases accord- 
ing to the Hagedom hypothesis [16, 17, 18, 19], such that their net effect turns out to 
be important. Indeed, hadronic resonances have been included in numerous applica- 
tions of the statistical models used in the studies of the ratios of hadron multiplicities 
[14, 20, 21, 22, 23] and the effects connected with their decays are essential for the 
successful description of the data. We note that only a quarter of the observed pions at 
RHIC comes from the "primordial" pions present at freeze-out, and the remaining three 
quarters are produced from the decays of resonances. 

In our approach, the same number of the resonances is included in the calculation of 
the hadron ratios and in the calculation of the hadron spectra. In this way, our theoretical 
spectra have always the correct relative normalization. Moreover, we have worked out 
semi-analytic formulas for the treatment of the resonance decays [4]. This allows us to 
sum up exactly many small contributions, especially, those appearing in the sequential 
decays. All two- and three-body decays are taken into account with the branching ratios 
taken from the tables. In the case of the three-body decays, the matrix elements are 
approximated by a constant, hence only the phase-space effect is included. 



Expansion 

In order to calculate the spectra we need to specify the expansion of the matter at 
freeze-out; clearly, the Doppler effect due to flow modifies the spectra and must be 
properly included. Our choice of the freeze-out hypersurface and the four-velocity at 
freeze-out has been made in the spirit of Refs. [11, 24, 25, 26, 27, 28, 29, 30] and is 
defined by the two conditions: 



t^ — r^ — rj — rj = const. , (2) 
and 

u>'=^ = Ui/JlX'A. (3) 
z z \ t t t J 

The constant in Eq. (2) will be later denoted simply by T. In order to make the transverse 
size, 

P = V^l+^, (4) 

finite, we impose the condition 

P < Pmax- (5) 

We note that the four-velocity (3) defining the hydrodynamic flow at freeze-out is pro- 
portional to the coordinate (Hubble-like expansion). This form of the flow and the fact 
the the coordinates t and are not limited and appear in the boost-invariant combination 
in (2) mean that our model is boost-invariant. In practical calculations it is convenient to 
introduce the following parameterization [28]: 

t = Tcosha|| coshcc^, r, = Tsinha|| cosh oCj^, 
rx = TsinhoCj^cos0, = T sinh ct^^ sin , (6) 

where tty is the rapidity of the fluid element, Vj- = r^/t = tanha||, and describes the 
transverse size, p = Tsinha^. The transverse velocity is Vp = tanhaj^/cosha||. The 
element of the hypersurface is defined as 

dx" dx^ dx^ 

where x^ = t,x^ = r^, x^ — ry, x^ = r^, and e^^^jg^ is the Levi-Civita tensor. A straight- 
forward calculation yields 

d'L^{x) = u'^{x) smh{a_Jcosh{a_Jda_^da^^d(l). (8) 

Equation (8) shows that the four- vectors dL^ and are parallel. In this case the spectra 
may be obtained from the expression analogous to the Cooper-Frye [31, 32] formula 



J^^ = jp^di:,f{p.u), 



(9) 



but with the distribution / which has collected the products of resonance decays (for 
details see [3]). With parameterization (6) we can rewrite Eq. (9) in the form 

dN o /■+°° /•Pmax/T _ r'iTl 



Inm^dm^dy 



= T dUu sinha^ J (sinha^j / d^p-uf{p-u), 

J —oa II JQ Jo 

(10) 



where 

p-u — mj^cosha||Coshaj^ — jc»^cos(^ sinhoCj^. (11) 

One can notice that the spectrum (10) is, as expected from the assumed boost invariance, 
independent of the rapidity y. 



RATIOS OF HADRON ABUNDANCES 

In the case of the boost-invariant systems, the ratios of hadron multiplicities at midra- 
pidity, dN/dy\y^Q, are simply related to the ratios of the local densities, rif, since 



dNjdy 



dNj/dy 



y=0 J J 



The first part of this equality follows from the boost invariance, whereas the second 
part is a consequence of the factorization of the volume of the system (this point is 
discussed in detail in Ref. [10]). Since the firecylinder formed at RHIC is approximately 
boost-invariant (at least within one unit of rapidity at y = 0, which is sufficient for 
our considerations which concentrate only on the central region), the ratios at zero 
rapidity measured for various particles may be used to fit the thermal parameters of 
the model. This is an important observation indicating that the ratios are not sensitive 
to the particular form of expansion. Consequently, the parameters of our model can be 
fixed in two steps: with help of the ratios we first fix the thermodynamic parameters, and 
later with help of the spectra we fix the two expansion parameters, T and Pmax- 
The density of the ith hadron species is calculated from the ideal-gas expression 



gJd'p/Hp), 



f^'^(P^ = 7^i^^P[iEiiP)-Mi-^sSi-^iIi)/T]±l) \ (13) 



where is the spin degeneracy, B-, S^, and 1^ denote the baryon number, strangeness, 
and the third component of isospin, and = ■\J p^ + mf. The quantities jUg,jU^ and 

llj are the chemical potentials which enforce the corresponding conservation laws. We 
recall that Eq. (13) is used to calculate the "primordial" densities of stable hadrons as 
well as of all resonances at freeze-out, which decay later on. 

Initially, the temperature, T, and the baryon chemical potential, fig, were fitted with 
the method to the 9 preliminary experimental ratios of the hadron yields (for the 
list of the ratios used in this fit and for more details of our analysis see Refs. [4, 5]). 
The fig and jlj were determined with the conditions that the initial strangeness of the 
system is zero, and the ratio of the baryon number to the electric charge is the same as 
in the colliding nuclei. This procedure gave usT = 165 ± 7 MeV and jUg = 41 ± 5 MeV. 
In Table I we present the result of the fit with a wider set of the up-to-date available 
hadronic ratios. It yields T = 168 ± 5 MeV and jU^ = 41 ± 4 MeV. It is interesting to 



TABLE 1. Optimal thermal parameters and the ratios of hadron 
multiplicities at zero rapidity. 



Model Experiment 



Fitted thermal parameters 

T [MeV] 168±5 

[MeV] 41 ±4 

Us [MeV] 10 

Hi [MeV] -1 

X^/n 0.6 

Theoretical ratios and the data 

1.02 1.00 ±0.02 [35], 0.99±0.02[36] 
0.09 0.08 ±0.01 [37] 

0.92 ±0.03 [38], 0.93 ±0.07 [39] 
^•^ 0.91 ±0.09 [35], 0.92 ±0.06 [36] 
0.16 0.15 ±0.02 [40] 

0.047 0.042 ±0.011 [13] 

0.042 0.039 ±0.011 [13] 

0.61 ±0.07 [37], 0.64 ±0.08 [39] 
0.60 ±0.07 [35], 0.61 ±0.06 [36] 
0.69 0.71 ±0.04 [38] 

0.77 0.83 ±0.06 [38] 

0.020 0.021 ±0.001 [41] 

0.15 0.13 ±0.03 [41] 

0.47 0.49 ±0.03 [42, 43] 

0.0011 0.0012 ±0.005 [44] 

0.0072 0.0088 ±0.0020 [45] 

0.86 0.95 ±0.16 [38] 



observe tliat tlie newly released data are fully consistent with the thermal picture and 
only a small change of the thermodynamic parameters follows when the new set of the 
ratios is used. 

A characteristic feature of our fit is that the optimal temperature is consistent with the 
value for the deconfinement phase transition obtained from the QCD lattice simulations: 
Tc = 154 ± 8 IVIeV for three massless flavors and Tc— 173 ± 8 iVleV for two massless fla- 
vors [33]). This type of the behavior has been also found in other statistical calculations 
[14, 34] and is interpreted as the argument for a scenario in which the hadronic ratios 
are fixed just in the hadronization process. 



TRANSVERSE-IVIASS SPECTRA 

Having fixed the two independent thermodynamic parameters of the model, we can 
use Eq. (10) to fit the transverse-mass spectra and to fix the two remaining geometric 
parameters of the model. This method has been initially applied in the study of the 



Tt /7r+ 
p/n- 

K-/K+ 

K- 1%- 

K*/h- 

P/P 

A/A 
Z/Z 

(^/K- 
A/p 
Qr/h- 
Z~ /n~ 




0.5 1 1.5 2 2.5 

m - m [GeV] 



FIGURE 1. The transverse-mass spectra at midrapidity. The data from STAR are denoted by asterisks, 
other symbols are used to denote the PHENIX data. All spectra are for most central collisions 



spectra of pions, kaons and protons, which gave us, for the most central events, the 
following values of the size parameters: 

T = 7.66 fm, pmax = 6.69 fm. (14) 



In Fig. 1 we show our results for all up-to-now available spectra at ^/s^ =130 GeV 
for the most-central collisions. In the upper part of Fig. 1 we show the spectra of 
pions, kaons, and antiprotons (used earlier to determine the geometric parameters) and 
the predicted spectra of the ^ and ^*(892)'^ mesons. The predicted spectrum of the 
mesons agrees very well with the measurement [41], with model curve crossing 
five out of the nine data points. The ^ meson production is of a particular interest in 
relativistic heavy-ion collisions, since its spectrum reflects the initial temperature of 
the hadronic system. This is because its interaction with the hadronic environment is 
negligible. Moreover, it does not receive any contribution from resonance decays. Thus, 
the agreement of the model and the data supports the idea of one universal freeze-out. 
The upper part of Fig. 1 also shows the averaged spectrum of K* resonances, with 
the data from Ref. [13]. Once again we observe a very good agreement between the 
model curve and the experimental points. As already mentioned above, the successful 
description of both the yield and the spectrum of ^*(892)^ mesons supports the concept 
of the thermal description of hadron production at RHIC, and brings evidence for a very 
small interval between chemical and thermal freeze-outs. 

In the lower part of Fig. 1 we show the predictions of the model for the spectra of 
hyperons. Again, in view of the fact that no extra parameters are introduced here and no 
refitting has been performed, the agreement is impressive. We note that the preliminary 
[46] data for the S's used in the figure were subsequently updated, which resulted in 
the reduction of the data by about a factor of 2. This correction makes our agreement 
with the data even better. The data accumulated at lower energies at SPS showed that the 
slope of the Q. hyperon was much steeper than for other particles [47]. On the contrary, 
in the case of RHIC the model predictions for the Q. are as good as for the other hadrons. 
Since the Q. contains three strange quarks, it is most sensitive for modifications of the 
simple thermal model used here, e.g. the use of canonical instead of the grand-canonical 
ensemble. The agreement of Fig. 1 does not support the need for inclusion of these 
effects. 



SUMMARY 

The success of our model in reproducing the hadron ratios and the transverse-mass 
spectra indicates that the particles are indeed produced thermally. The model has only 
four parameters and describes the data with surprising accuracy. Since our approach uses 
hadronic degrees of freedom and starts at freeze-out, important theoretical questions 
concerning the earlier stages of the evolution cannot be addressed in this framework. 
We think, however, that our results constrain any more microscopic descriptions of 
the evolution of the matter formed in ultra-relativistic heavy-ion collisions. Further 
applications of the model aiming at the description of the elliptic flow and the HBT 
pion radii will be presented in the next contribution to these Proceedings [10]. 
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